SMALE'S FUNDAMENTAL THEOREM OF ALGEBRA 

RECONSIDERED 



DIEGO ARMENTANO AND MICHAEL SHUB 

Abstract. In his 1981 Fundamental Theorem of Algebra paper 
Steve Smale initiated the complexity theory of finding a solution of 
polynomial equations of one complex variable by a variant of New- 
ton's method. In this paper we reconsider his algorithm in the light 
of work done in the intervening years. Smale's upper bound esti- 
mate was infinite average cost. Our's is polynomial in the Bczout 
number and the dimension of the input. Hence polynomial for any 
range of dimensions where the Bczout number is polynomial in the 
input size. In particular not just for the case that Smale considered 
but for a range of dimensions as considered by Biirgisser-Cucker 
where the max of the degrees is greater than or equal to n 1+e for 
some fixed e. It is possible that Smale's algorithm is polynomial 
cost in all dimensions and our main theorem raises some problems 
that might lead to a proof of such a theorem. 



1. Introduction and Main Result 

1.1. Introduction and Preliminaries. In his 1981 paper [20] Steve 
Smale initiated the complexity theory of finding a solution of poly- 
nomial equations of one complex variable by a variant of Newton's 
method. More specifically he considered the affme space of monic 
polynomials of degree d, 
d 

f{z) = s ^a i z\ a d = l and a, G C, (i = 0, . . . , d — 1). 

i=0 

He identified P^ with C d , with coordinates (do, • • • , dd-i 

) G C d . In P d 

he considered the polydisk 

Qi = {/GP d : H < 1, i = 0,...,d-l} 
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to have finite volume and he obtained a probability space by normal- 
izing the volume equal 1. The algorithm he analyzed is given by: let 
< h < 1 and let z = 0. Inductively define z n = T h (z n _i) where T h is 
the modified Newton's method for / given by T^z) = z — hjf^j. 
His eponymous main theorem was: 

Main Theorem: There is a universal polynomial S(d, 1/ fi) 
and a function h = h(d, fi) such that for degree d and \i, 
< fi < 1, the following is true with probability 1 — fi. 
Let xq = 0. Then x n = Th(x n _i) is defined for all n > 
andx s is an approximate zero for f where s = S(d, 

In [20], that x s is an approximate zero meant that there is an x* 
such that f(x*) = 0, x n —> x* and ^§*yp < |, for j > s, where 

Xk+i = Xk — ■ That is, Xk+i is defined by the usual Newton's 
method for /. Smale mentions that the polynomial S may be taken 
to be 10 °(^+ 2 ) _ xhe notion of approximate zero was changed in later 
papers (see Blum et al. [H] for the later version or Section The 
new version incorporates immediate quadratic convergence of Newton's 
method on an approximate zero. In the remainder of this paper an 
approximate zero refers to the new version. 

Note that -jj is not finitely integrable, so Smale's initial algorithm 
was not proven to be finite average time or cost when the upper bound is 
averaged over the polydisk Qi (see Blum et al. [U page 208, Proposition 
2])- 

A tremendous amount of work has been done in the last 30 years 
following on Smale's initial contribution, much too much to survey 
here. Let us mention a few of the main changes. In one variable a lot 
of work has been done concerning the choice of good starting point zq 
for Smale's algorithm other than zero. See chapters 8 and 9 of Blum 
et al. (8] and references in the commentary on chapter 9. The latest 
work in this direction is Kim-Martens-Sutherland [T3] . 

Smale's algorithm may be given the following interpretation. For 
z e C, consider f t = f - (1 - t)f(z ), for < t < 1. f t is a poly- 
nomial of the same degree as /, z is a zero of f and fi = f. So, we 
analytically continue Zq to z t a zero of f t . For t — 1 we arrive at a zero 
of /. Newton's method is then used to produce a discrete numerical 
approximation to the path (f t ,z t ). 

If we view / as a mapping from C to C, then the curve Zt is the branch 
of the inverse image of the line segment L = {tf(zo) '■ < t < 1}, 
containing z . (See Figure^ . 
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Figure 1. 



Here are several of the changes made in the intervening years. Rene- 
gar [13] considered systems of n-complex polynomials in n-variables, 
without the restriction to be monic. Given a degree d, we let stands 
for the vector space of degree d polynomials in n complex variables 

7 d = {/ : f(z) = 

\\a\\<d 

where a = (cki, . . . , a n ) 6 N" is a multi-index, ||a|| = Ylt=i z °" = 
Z V~ ' ' ' z n n i a a ^ C. We have suppressed the n for ease of notation. It 
should be understood from the context. 

For (d) = (di, . . . , d n ), let ? {d) = ? dl x • • ■ x ? dn so / = (f u . . . , /„) G 
J^d) is a system of n polynomial equations in n complex variables and 
fi has degree di. 

As Smale's, Renegar's results were not finite average cost or time. 
In a series of papers Shub and Smale [TB], [T7], [IB], [H], made some 
further changes and achieved enough results for Smale 17th problem 
to emerge a reasonable if challenging research goal. Let us recall the 
17th problem from Smale [21] : 

Problem 17: Solving Polynomial Equations. 

Can a zero of n- complex polynomial equations in n-unknowns 
be found approximately, on the average, in polynomial 
time with a uniform algorithm? 

In place of 7(d) h is natural to consider W(d) — Hdi x • • • x %d n i 
where r H.d i is the vector space of homogeneous polynomials of degree di 
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in n + 1 complex variables. The map 

hi '■ — > Hdf, irfi(/)(^0) • • • j z n) = Zq /(—,...,— ) , 

is an isomorphism and i : 7^) H(d) f° r i = (idi> • • • , id„) is an isomor- 
phism. 

For / G and A G C, 

/(AC) = A (A*) /(C), 

where A(oj) means the diagonal matrix whose z'-th diagonal entry is <ij. 
Thus the zeros of / G H(d) are now complex lines so may be considered 
as points in projective space P(C n+1 ). 
The affine chart 

j : C n P( C - +1 ), j(d, . . . , Cn) = C(l : Ci : • • • : Cn), 

maps the zeros of / G 7^) to zeros of i(/) G %{d)- In addition i(/) may 
have zeros at infinity i.e. zeros with Co = 0. 

From now on we consider %{d) and P(C n+1 ). On T-L^ we put a unitar- 
ily invariant Hermitian structure which we first encountered in the book 
[22] by Hermann Weyl and which is sometimes called Weyl, Bombieri- 
Weyl or Kostlan Hermitian structure depending on the applications 
considered. 

For a = (oto, ■ ■ ■ , a n ) G N n+1 , ||a|| = di, and the monomial z a = 
z o° ' ' ' z n n the Weyl Hermitian structure makes (z a , z@) = 0, for a ^ (3 
and 

'di\~ ( di\ 



(z a ,z a ) 

\a J \ocqI ■ ■ ■ a r 
On H{d) we P u t the product structure 



(f,9) = ^2ifi,9i)- 



i=i 



On C we put the usual Hermitian structure 

n 

fay) = ^2xkVk- 

Given a complex vector space V with Hermitian structure and a 
vector ^ v G V, we let v L be the Hermitian complement of v, 

v 1 ' = {w G V : (v,w) = 0}. 

v 1 - is a model for the tangent space, T V F(V), of the projective space 
P(V) at the equivalence class of v (which we also denote by v). 
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T V F(V) inherits an Hermitian structure from (•, •) by the formula 

{w 1 ,w 3 )v = —, r— , 

where W\, W2 G v 1 - represent the tangent vectors at T V ¥(V). 

This Hermitian structure which is well defined is called the Fubini- 
Study Hermitian structure. 

The group of unitary transformations U(n+1) acts on Hid) and C n+1 
by / / o U- 1 and C H> U( for U G W(n + 1). 

This unitary action preserves the Hermitian structure on Hu) and 
C n+1 , see Blum et al. 0. That is, for U G W(n + 1), 

(foU-\goU- 1 ) = (f,g) for /, <? G H(d)\ 

(U(,U(') = {(,C) for C,C eC" +1 . 

The zeros of Xf and / for ^ A 6 C are the same, and we may 
consider the space P(ft(d)). Now the space of problem instances is 
compact and the space P(C n+1 ) is compact as well. P {H(d)) has a 
unitarily invariant Hermitian structure which gives rise to a volume 

N— 1 

form of finite volume ytnJj where N = dimH^)- 
The average of a function <\> : P {H(d)) — > K. is 

Ew« vU) = -4 rr / Mf)df = I <j>(f)df. 

If is induced by a homogeneous function : Hid) — > R of degree 
zero, that is, <fi(\f) = <f>(f), A G C — {0}, then we may also compute 
this average with respect to the Gaussian measure on (ttfd), (■, ■)), that 
is, 

(1) S,^)^-/ 0(/)e-^ 2 / 2 rf/. 

And it is this approach via the Gaussians above defined on T-i(d) and 
the Fubini-Study Hermitian structure and volume form on P(C n+1 ) 
that we take in this paper. The quantities we define on %(d) are homo- 
geneous of degree zero, thus are defined on P [H(d)) and benefit from 
the compactness of this space and of P(C n+1 ). While averages over 
systems of equations may be carried out in the vector space H(d)- 

The solution variety 

V = {(/,*) G {U [d) - {0}) x P(C" +1 ) : f{x) = 0}, 
is a central object of study. 
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V is equipped with two projections 

V 





(d) 



P(C" +1 ) 



The solution variety V also has a projective version, namely, 
V P = {(/, x) e P (n (d) ) x P(C" +1 ) : f(x) = 0}. 

1.2. Homotopy Methods. Homotopy methods for the solution of 
a system / e proceed as follows. Choose (gf, () £ V a known 
pair. Connect g to / by a C 1 curve / f in "H^), < t < 1, such that 
fo — 9i fi — /) an d continue Co = C to Ct such that ft((t) = 0, so that 
/i(Ci) = 0- The critical values of the projection of V on %(d) — {0} 
are an algebraic subvariety, E, of %(d) — {0} of complex codimension 
1, called the discriminant variety. By the transversality theorem (see 
Abraham- Robbin [2]) a generic set of C 1 curves f t do not intersect E. If 
a curve is in this generic set and /(Co) — 0, then by the implicit function 
theorem we may continue Co to Ct, < t < 1, such that ft((t) = 0. See 
Smale [20] for this type of argument. Indeed almost all "straight line" 
paths in Hu) do not intersect E, again by a transversality argument, 
so if Co is a nondegenerate zero of g then for almost all /, Co m ay be 
continued to a zero of / along the curve f t — (1 — t)g + tf. We do not 
use this generality in this paper so we leave the above assertions as a 
sketch. In Proposition[7\we prove a precise version of the fact that the 
homotopies we consider in this paper may be almost always continued. 

Now homotopy methods numerically approximate the path (/t, Ct)- 
One way to accomplish the approximation is via (projective) Newton's 
methods. Given an approximation x% to Ct define 

Xt+At = N ft+At (x t ), 

where 

N f (x)=x-(Df(x)U)- 1 f(x). 

That Xt is an approximate zero of ft associated with the zero Ct means 
that the sequence of Newton iterative Nf t (xt) converges immediately 
quadratically to Ct- 

The main result of Shub [15] is that At may be chosen so that to = 0, 
tk = t k _i + At k , x tk is an approximate zero of f tk with associated zero 
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(t k , and t K = 1 for 

(2) K = K(J, g,C)<C D 3/2 f n(f t , Ct) || (/*, Ct) || c/ t A ) 
Here C is a universal constant, D = maxdj, 

M/,C) = 11/11 •||W(OI^)" 1 ^(IICII*- 1 >/*)ll 

is the condition number of / at (, and 

ll(/U*)ll(/ t , Ct ) = (ll/*ll/, + IIC i |lc t ) 1/2 

is the norm of the tangent vector to the projected curve in (ft,Ct) m 
V P C P (n {d) ) x P(C n+1 ). The choice of At k is made explicit in Dedieu- 
Malajovich-Shub [H] and Beltran [I]. 

In Vp, HCtllcf — A*(/t)Ct) ||/*||/t) so the estimates (j2J) may be bounded 
from above by 

(3) K(f,g,C)<CD 3 / 2 [\(f t Xt) 2 \\ft\\f t dt, 

Jo 

for a perhaps different universal constant C. 

Finally in the case of straight line homotopy ||/< ||/ t = - 777^2 ' 

\\jt\\ 

where 9 is the angle between /o and /1. So ([3]) may be rewritten as 

(4) K(f,g,C) < CD 3 / 2 Edn(fl) ||/o|| ||/i|| f dt, 

see Biirgisser-Cucker [9j, where (jl]) is a principal part of the analysis 
and where the increments A tk which exhibit the right hand side of @ 
as an upper bound are also made explicit. 

Much attention has been devoted to studying the right hand of (jl]), 
for a good starting point (g,C)- 

In Beltran-Pardo [5], an affirmative probabilistic solution to Smale's 
17th problem is proven. They prove that a random point (g, £) is good 
in the sense that with random fixed starting point (g, () = (/o, Co) 
the average value of the right hand side of (jl]) is bounded by 0(nN). 
Moreover, Beltran and Pardo show how to pick a random starting point 
starting from a random n x (n + 1) matrix. 

In [9], Biirgisser-Cucker exhibit a deterministic algorithm for Smale's 
17th problem which is polynomial average cost, except for a narrow 
range of dimensions. More precisely: 

There is a deterministic real number algorithm that on 
input f G Hid) computes an approximate zero of f in 
average time N°^ oglogN \ where N = dim H(d) measures 
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the size of the input f . Moreover, if we restrict data to 
polynomials satisfying 

D < n~^, or D > n 1+£ , 

for some fixed e > 0, then the average time of the algo- 
rithm is polynomial in the input size N . 

So Smale's 17th problem in its deterministic form remains open for a 
narrow range of degrees and variables. 

The case D < fjW is dealt with by Biirgisser-Cucker by constructing 
a good starting point for a homotopy method while the case D > n 1+£ 
is dealt with differently. Our Theorem [3] shows that we may use the 
homotopy method suggested by Smale's algorithm and described in 
the next section in this range of dimensions and conclude a polynomial 
result as well. 

1.3. Smale's algorithm reconsidered. Smale's 1981 algorithm de- 
pends on /(0), so there is no fixed starting point for all systems. Given 
( G P(C n+1 ) we define for / G Htd) the straight line segment f t G Hu), 
< t < 1, by 

/« = /-(l-i)A(|||)/(C). 

So /o(C) — and f\ = f. Therefore we may apply homotopy methods 
to this line segment. (Here A ^'^'. ^ must be understood as a ma- 
trix of functions and f(() as a vector of constants which multiply the 
functions according to matrix vector multiplication.) 

Note that if we restrict / to the affine chart ( + £ then 

/,(*) = /(*)-(i-f)/(C), 

and if we take ( = (1, 0, . . . , 0) and n = 1 we recover Smale's algorithm. 

There is no reason to single out ( = (1, 0, . . . , 0). Since the unitary 
group acts by isometries on P ("H(d)), P(C n+1 ), V and Vp, and preserves 
/i and is transitive on P(C n+1 ), all the points ( yield algorithms with 
the same average cost. 

Note that if we let 

V c = {/ G H (d) : /(C) = 0}, 

then 

/„ = /-A(|gl)/(C), 

is the orthogonal projection n^(/) of / on V<^. This follows from the 
reproducing kernel property of the Weyl Hermitian product on 
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namely, 

(5) <0,<-,C>*> = <7(O, 

for all g G (i — 1, . . . , n). In particular || (-, (} di \\ = \\(\\ di - 
Then, 

||/-nc(/)ll = l|A(||CI|-*)/(C)ll, 

while 

||n c (/)|| = (||/|| 2 -||A(||ci|- d 0/(C)ll 2 ) 1/2 . 

Let $ : n (d) x P(C" +1 ) x [0, 1] ->■ V be the map given by 

(6) *(/,C,*) = (/*,&)» 
where 

/ t = (i-t)n f (/) + */, 

that is, 

/* = /-(l-t)A(|^)/(C), 
and is the homotopy continuation of ( along the path f t . 

Proposition 1. For almost every f G %{d)> the set of C £ P(C n+1 ) 
snc/i £/iat $ defined for all t G [0, 1] /ias full measure. Moreover, for 
every ( G P(C ra+1 ), the set of f G %r<2) suc/i that $ zs defined for all 
t G [0, 1] has full measure. 

Remark: In fact, the proof also shows that the complement of the 
set (/, () such that $ is defined for all t G [0, 1] is a real algebraic set. 
The proof of Proposition U\ is in Section [B 

The norm of f t is given now by the formula 

m _ Il/oll \\h\\ MQ) _ PcCQJ ll/-n c (/)|| 
IIJ<II/ * ll/*ll 2 ll/*ll 2 

(ii/ii 2 -iiA(iicir d ')/(oii 2 ) 1/2 iiA(iicir*)/(oii 

WftW 2 

Let r(/,C) = ^(/,n c (/),C) and T € {f) = T(f,()- Then, the average 
cost of this algorithm satisfies 

Proposition 2. 

E W(d) (r c )=E W(d)XP(c „ +1) (T)<(I), 

where 



(2 7 r)^ w /(p(c™+ 1 )) y /6Ww) y ceP(Cn+1) y te[0)1] h/j 2 
-iiA(iicir d< )/(c)ir) 1/2 iiA(nci|-*)/(c)iie- ||/||2/2 d/dc^. 
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As we have mentioned above it is easy to see by unitary invariance of 
all the quantities involved that the average E-^. (7c;), on H^, is inde- 
pendent of C and equal to E W(d)XF ( C n+i)(T). This argument proves the 
first equality of this proposition. The inequality follows immediately 
from the definition of T(f, C). 

What is gained by letting £ vary and dividing by vol(P(C n+1 )) is that 
a new way to see the integral leads to new theorems and interesting 
questions. 

Suppose 77 is a non-degenerate zero of h G T-Lu)- We define the 
basin of 77, B(h,r)), as those £ G P(C rt+1 ) such that the zero £ of 

h — A ( ^'^i ) MC) continues to 77 for the homotopy h t . From the proof 
of Proposition[J\we observe that the basins are open sets. 

Let (I) be the expression defined in Proposition^ Then, the main 
result of this paper is: 

Theorem 1 (Main Theorem). 



Hd) r)/h(rf)=0 



where 



6(M) = 7^Im / **(CK. 



7J0Z(P(O +1 )) 

e fc (C) = (||/i|| a -||A(||ci|-*)MC)ll 2 ) 1/2 - 

• iiA(iicir^)Mc)ii-x n (iiA(iicir^(0H72), 

and I„(a) = eCi-*" 2 )^- 2 "- 1 ^. 

From Proposition [7] we obtain that the function O, defined in the 
statement of TheoremUl is defined for almost every pair (h,rj) G V. 

Summing Q(h, 77) over the roots of h we let Q(h) = Y2n/h(rf)=o ®(^' v)j 
and for almost all h we have that 

voi(P(c n+i )) y CeP( c«+i) 

That is, Q(h) = \\e h \\ L i. 

More generally, for p > 1, consider the L p -norm of &h 
( 8 ) II^IIlp = / ^(C) p dC 

voi(P(c n+i )) y CeP( c«+i) 

The next theorem shows that the average of H^H^p over H^, for all 
p > 1, is polynomial in AT. 
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Theorem 2. 



„ 2P T(N-n + p/2) T(n + p/2) 



^H ld) (\\9h\\ L p) — T^T7 \ FT - , 

w p I (iV — n) 1 (n) 

(The equality holds for p = 1). 
Theorem 3. 

(I) < IZCD^Vn^N^ 2 . 

That is, (I) is polynomial in the Bezout number and the input size, 
N, and polynomial in the input size alone for any range of dimensions 
where the Bezout number D is polynomial in N. 

Since our method of proof of Theorem [3] relies on Theorem where 
the basins are not taking into account, it is possible that Smale's algo- 
rithm is polynomial cost in all dimensions. 

Understanding the basins better might lead to a proof of such a 
theorem. The integral 

V /AM) M y 2 < e(n + l) 
2^ 117,112 6 an^ v, 



{2k] 



N 



where T) — d± ■ ■ ■ d n is the Bezout number (see Biirgisser-Cucker [9]). 
So the question is how does the factor 0(h,r)) affect the integral. □ 

From Theorem® the expected value of 0(h) = WO^Wl 1 is controlled, 
then, if the integral on the D basins are reasonably balanced, the factor 
of D in Theorem and the integral above may cancel. 

Remark: The proof of Theorem [J\ involves complicated formulas 
which exhibit enormous cancellations. We do not have a good ex- 
planation for these cancellations. 

At the end of the paper we present some numerical experiments 
with n = 1 and d = 7 which were done by Carlos Beltran on the 
Altamira super computer at the Universidad de Cantabria (partially 
supported by MTM2010-16051 Spanish Ministry of Science and Inno- 
vation MICINN). We thank Carlos and the Universidad de Cantabria. 
We also thanks Gregorio Malajovich for many useful discussions and 



1 In an earlier version of this paper we asked: 
(d) Evaluate or estimate 

1 



cep ( o +1 ) IIAdiciMOMOII 2 ' 1 - 1 



e iiA(iicir d *)Mc)ii 2 /2 d( - 



It is easy to see, as in the proof of Theorem 2, that the expected value of this 
integral over H(d) is infinite. In Fernandez-Pardo [T2] the authors consider the 
more meaningful average over the unit sphere of 'H(d) and get a precise formula for 
it. Our initial goal in asking question (d) was to get an upper-bound estimate of 
the integral we now evaluate in Theorem 2. 
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Santiago Laplagne for having done some more experiments. It would 
be interesting to see more experimental data. The proof of Theorem 
U\ is in Section and the proofs of Theorem [H and Theorem [3] are in 
Section [^J 

2. Proof of Proposition [TJ 

For the proof of Proposition [J\ we need a technical lemma. 

Lemma 1. Let E be a vector bundle over B, F be finite dimensional 
vector space, and consider the trivial vector bundle F x B. Let ip : 
F x B E be a bundle map, covering the identity in B, which is a 
fiberwise surjective linear map. Then, if is a surjective submersion. 

The proof is left to the reader. 

Recall that $ : U {d) x P(C n+1 ) x [0, 1] V is the map given by 

mC,t) = (ft,Ct), 

where 

/* = /-(l-*)A(^)/(C), 

and is the homotopy continuation of ( along the path f t . 

This map is defined at (/, (,t) provided that rank(D/ s (C s )|^±) = n, 
for allj^e [0,t]. 

Let K be the vector bundle over C n+1 — {0} with fiber K z = L(z ± , C n ), 
z G C" +1 — {0}, where L(z ± , C n ) is the space of linear transformations 
from z 1 - to C n . For k — 0, . . . , n, let Kk be the sub-bundle of rank k 
linear transformations. Note that Kk has (n — k) 2 complex codimension 
(c.f. [3]). These sub- bundles define a stratification of the bundle K. 

Lemma 2. Let fi (0) be the set of pairs (f, () G H{$ x P(C n+1 ) such 
that rank{Df Q {C)\ c x) < n, where (/ ,C) = *(/,C,0). Then is a 
stratified set of smooth manifolds of complex codimension (n — k) 2 , for 
k = 0, . . . , n — 1. 

Proof. Let be the inverse image of VL^ under the canonical pro- 
jection H (d) x (C n+1 - {0}) -> x P(C n+1 ). 

Let y> : H (d) x (C n+1 -{0}) -> K be the map </?(/, C) = Df(()\ c ±. For 
each /c = 0, ... let 6< 0) = ^(iQ. Since D/ (C)I^ = Af(C)l^, 

then = U™- 1 ^ 0) . 

Claim: if is transversal to Kk for = 0, . . . , n — 1: 
Note that <p(f, ■) : C n+1 — {0} — > K is a section of the vector bundle 
iY for each / G "H(d). Moreover, for each ( G C n+1 — {0}, the linear 
map ip(-, C) : "H(d) — >■ -^f is a surjective linear map. This fact follows by 
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construction: given L e~K c = L(C X , C n ), let L G L(C n+1 ,C n ) be any 
linear extension of L to C n+1 . Then, the system / = A( ^'^ d '._ 1 )L(-) G 

H(d) satisfy Df(Q\^± = L. Then, the claim follows from Lemma QJ 

Since tp is transversal, we conclude that the inverse image of a strat- 
ification is a stratification of the same codimension (c.f. [3]). That is, 
is a stratification of complex submanifolds of complex codimension 
{n — k) 2 , for k — 0, . . . , n — 1. 

Moreover, since each strata fij^. contains the fiber of the canonical 
projection U( d ) x (C n+1 - {0}) U {d) x P(C n+1 ), then, its image, 

ilk , is a smooth manifold of codimension (n — /c) 2 , and the lemma 
follows. □ 

One can define the homotopy continuation of the pair (/, () G Wte) x 
P(C n+1 ) for all t G [0,1] whenever (f,Q £ and lies outside the 
subset of pairs such that there exist (w, t) G P(C™ +1 ) x (0, 1] satisfying 
the following equations: 



Note that, since f t is homogeneous, rank(D f t (w)\ w ±) is well defined on 
w G P(C n+1 ). 

Differentiating f t we get 



Therefore, taking s = 1 — t, we conclude that one can define the homo- 
topy continuation of the pair (/, C) G H($ x P(C n+1 ) for all t G [0, 1] 
whenever (/, Q and lies outside the subset of pairs such that 

there exist (w, s) G P(C n+1 ) x [0, 1) satisfying, for some k — 0, . . . , n— 1, 
the following equations: 



rank ( [A«£ C>*) ■ Df(w) - s A(di(w, C)*"^-, C» /(C)] D = k. 



Let E' C V be the set of critical points of the projection m : V — > 
%{d)- Recall that E = 7Ti(E') C H(d) is the discriminant variety. If 
(/, w) G E' then w is a degenerate root of /, that is, r&nk(D f (w)\ w ±) < 
n (c.f. Blum et al. 0). 

Note that if / G E then every ( G P(C n+1 ) satisfies equations © 
and (fTUj) for s = and w G P(C n+1 ) a degenerate root of /. Hence, 





(9) 



a«c,C>*)/M - s A(( w ,c>*)/(0 = o, 



(10) 
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it is natural to remove the discriminant variety £ and the case s = 
from this discussion. 

Lemma 3. Let A C {H {d) - E) x P(C n+1 ) x P(C n+1 ) x (0, 1) be the 
set of tuples (f, (, w, s) such that equations (Tj|) and [Td\) holds for some 
k = 0, . . . , n — 1. Then, A is stratified set of smooth manifolds of real 
codimension 2{n + (n — k) 2 ) for k — 0, . . . ,n — 1. 

Proof. Similar to the preceding proof, for each k = 0, . . . , n — 1, we 
consider the set A fc C {U [d) - £) x (C n+1 - {0}) x (C n+1 - {0}) x (0, 1) 
of tuples (/, C, w, s) such that equations ([9]) and ( TTDT) holds . 

Let (/, C) w, s) G A k for some k G {0, . . . , n — 1}. Since / ^ £ then 
(w,C) 7^ 0. Therefore from ([9]), equation ( |T0i) takes the form 

rank (((w, ()Df(w) - A(d<) /(«;) (•, 0) | tt x) = fc, 

for A; = 0, . . . , n — 1. 
Let 

F=(F 1 ,F 2 ) : (^ (d) -S)x(C n+1 -{0})x(C n+1 -{0})x(0,l) -»• C n xK 
be the map defined by 

F x (f, C, s) = A«C, C) dl )/H - s A« Wj C) di )/(C) G C n 
F 2 (/, C, s) = (w, ««;, C)5/(w) - A(dO /H (•, C» Ux) G if. 

Note that A fe = F _1 ({0} x TQ. 
Claim: F is transversal to {0} x K k . 

In fact, what we prove is that DF is surjective at any point (/, C, w, s) 
which F maps into {0} x K k , for any k = 0, ... ,n — 1, that is, any 
point in Afe. 

Recall that = {/ G 7/(d) : /(C) = 0}. Consider the orthogonal 
decomposition = V( © C^, where C{ = V^ 1 . 

Let (/, C, w,s) G A k . We first prove that DF x (f, C, w, s)|c { : — >■ 
C n is surjective. 

Note that the linear map £ : C n — >■ given by £(a) = A( ^'^ d '. )a, 

is an isomorphism, where £ _1 : — )■ C n is given by £ _1 (/) = /(C)- 
Then, under this identification, the restriction to of the derivative 
of Fx is the linear map given by 

DF 1 (f,(,w,s)\ c< = (1- s)A{{w,() di ), 

for all tuples (/, C, w, s). Moreover, since (/, C, w, s) G A fc , then (w, C) ^ 
and s ^ 1, hence DF 1 (f,(,w,s)\ r is onto. 
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Now we prove that DF 2 (f, C, w, s)\ B ^n+n is surjective onto the 

tangent space Tp 2 ^ jW>s ^K, at every (f,£,w,s) G A*. 

Note that the map F 2 (f, C, •, s) : C n+1 - {0} -»■ K is a section of the 
vector bundle K. Therefore, from LemmaUl it is enough to prove that 
F 2 (-, C, s) : H(d) — > -ft'w is a surjective linear map. 

Fix a tuple (/, £, w, s) G A fc , for some k = 0, . . . , n — 1. The unitary 
group W(n+1) acts by isometries on (% (d) -E) x (C n+1 -{0}) x (C n+1 - 

{0}) x (0, 1) by [/ ■ (/, C, «J, s) = (/ o U~\ U(C), U(w) 7 8), and leave A fc 
invariant. Therefore we may assume that w = cq. Write fi{z) = 
Y,\\a\\=di a <xZ a , (i = 1, ...,n). Then, the linear map F 2 (-,(,e , s) : 
H( d ) ->■ ^ eo is given by 

F 2 (f, C,eO,s) = ((CoO^i^) - ^ a (d 1 ,0,...,0)0))ij=l,.-,n, 

where u,- is the n-vector with the j-entry equal to 1 and the others 
entries equal to 0. 

In particular, since £ ^ 0, the restriction F 2 (-, C, Co, s) : Vf -> K eo is 
surjective, concluding the claim. 

Then, since F is transversal to {0} x K^, we conclude that A& = 
i 71_1 ({0} x Kf.) is a submanifold of real codimension 2(n + (n — fc) 2 ), 
for A; = 0, . . . , n — 1. 

To end the proof, we note that A& is transversal to the fiber of the 
canonical projection {U {d) - E) x (C n+1 - {0}) x (C n+1 - {0}) x (0, 1) -> 
(U (d) — E) x P(C n+1 ) x P(C n+1 ) x (0, 1). □ 

Let IT : H [d) x P(C n+1 ) x P(C rt+1 ) x (0, 1) U [d) x P(C n+1 ) be the 
canonical projection 

n(/,c,^, s ) = (/,C). 

Then, from Lemma d and Lemma the set of pairs (/, Q G "H(d) x 
P(C n+1 ) such that the homotopy is not defined for all t G [0, 1] is 
contained by the union 

Q {0) U n(A) U E x P(C n+1 ) C U (d) x P(C n+1 ). 

Remark: We could conclude the proof by Fubini's Theorem. But we 
give a different argument. See the remark at the end. 

Proof of Proposition^ For k = 0, . . . , n - 1, let Vt^ ] C H {d) x P(C n+1 ) 
be the subset given in the proof of Lemma [H and let t\\ : H( d ) x 
P(C n+1 ) — > T~L<d) be the projection in the first coordinate. From Sard's 
Lemma we get that almost every / G W(d) is a regular value of the 
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restriction 7Ti| n (o) : fl\ — > H(d), for each k — 0, . . . ,n — 1. Therefore, 
from Lemma we conclude that for almost every / E%(d) the subset 

^i| n (o," 1 (/) = 7 rr 1 (/)n^ 0) cP(C" +1 ), 

k 

is an empty set or a smooth submanifold of complex dimension n — 
{n — k) 2 , for k = 0, . . . , n — 1. Hence, for almost every / G "H(d), the set 
of £ G P(C n+1 ) such that $ is not defined at t — has measure zero. 

Similar to the preceding argument, for each k = 0, ... ,n — 1, let 
A fc C (H {d) - S) x P(C n+1 ) x P(C n+1 ) x P(C n+1 ) x (0, 1) be the set 
of tuples (/, C, w, s) such that equations (jUj) and (fTOj) hold, and let 
flf : ft (d) x P(C n+1 ) x P(C n+1 ) x P(C n+1 ) x (0, 1) H [d) be the 
projection in the first coordinate. Then by Sard's Lemma, almost every 
/ G H(d) is a regular value of the restriction flf\\ k '■ ^ k — > H^)- 
Therefore, from Lemma we conclude that for almost every / e H(d) 
the subset 

n f \ A ~\f) = nf(f) n A fc c p(c n+1 ) x p(C" +1 ) x (o, i), 

is an empty set or a smooth submanifold of real dimension 2n + 1 — 
2{n — k) 2 , for k = 0, ... ,n — 1. Then, projecting in the ^-space we 
obtain that for almost every / G %{d)i the set of ( G P(C n+1 ) such 
that $ is not defined at t G (0, 1) is a finite union of measure zero 
sets. Moreover, since £ C Hid) has measure zero, the proof of the first 
statement of the proposition follows. 

The second statement of Proposition [I] follows directly from proofs 
of the claims of Lemma [H and Lemma 0, and the subsequent analysis 
of dimensions. □ 

Remark: The proof of Proposition Q] follows immediately from Fu- 
bini's Theorem. But we say more because this discussion may be use- 
ful for the discussion of the basins. This proposition proves that the 
boundary of the basins are contained in this stratified set, the structure 
of which should be persistent by the isotopy theorem (c.f. j3]) on the 
connected components of the complement of the critical values of the 
projection. We don't know if there is more than one component. 

3. Proof of Theorem rjj 

Let us first state the notation in the forthcoming computations. 
Most of the maps are defined between Hermitian spaces, however they 
are real differentiable. Therefore, unless we mention the contrary, 
all derivatives are real derivatives. Moreover, if a map is defined on 
P(C n+1 ) then is natural to restrict its derivative at ( to the complex 
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tangent space T^P(C n+1 ). If L : E — > F is a linear map between finite 
dimensional Hermitian vector spaces, then its determinant, det(L), is 
the determinant of the linear map L : E — > Im(L), computed with 
respect to the associated canonical real structures, namely, the real 
part of the Hermitian product of E and the real part of the inherited 
Hermitian product on Im(L) C F. The adjoint operator L* : F — > E is 
also computed with respect to the associated canonical real structures. 

In general, if E is a set, Id# means the identity map defined on that 
set. 

Since the set of triples (/,C,*) e H(d) x P(C n+1 ) x [0,1] such that 
t = or t = 1 has measure zero, we may assume in the rest of this 
section that i 6 (0,1). 

Recall that $ : H {d) x P(C n+1 ) x [0, 1] -> V is the map given by 

*(/,C,*) = (/t,Ct), 

where 

/* = /-(i-*)A(|^)/(C), 

and ( t is the homotopy continuation of ( along the path f t . 

For each t G (0, 1), let $ t : H w x P(C n+1 ) -> V be the restriction 

**(■,-) = *(■>■.*)■ 

Recall that for each non-degenerate root r\ of h, B(h,rj) is the non- 
empty open set of those ( G P(C n+1 ) such that the zero ( of H$(h) 
continues to r] for the homotopy h t = (1 — t)H^(h) + th. 

Given h G U {d) and t G (0, 1), let H t : P(C n+1 ) U( d ) x P(C n+1 ), 

(11)^(0 = ^(0,0, and fc i (0 = fc+(^)A(M^)/i(0 ) 

for all £ G P(C n+1 ). (We have suppressed the h for case of notation.) 

Lemma 4. Let t G (0, 1), and let (h,rj) G V 6e a regular value of $ t . 
Then, the fiber Q t {h,r])~ l is given by 

^ l (h, V ) = H t (B(h, V )). 

Proof. For < t < 1, we have that (/, £ &7 1 (hi r l) provided that 

i) /i = / t = t/ + (l-t)nf(/); 

ii) the homotopy continuation of ( on the path {sh+(l— s)n^(/)} se r 0i i] 
is rj. 

Since n^(/i) = n^(/) we conclude that 

/ = ] (h - (1 - *)n c (/0) = + (^) (h - n c (h)), 

and C G £(M)- □ 
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Proposition 3. Let (f, () G H {d) x P(C n+1 ) such that $ t is defined 
and let (h,rj) = § t {f\0- Then the normal jacobian of $ t is given by 



NJ* t (f,0=t 



2n JaC H t (0 

NJ ni (h jV )' 



where Jac^ (() = \ det(DH t (())\ is the jacobian of the map H t defined 

m m\). 



The proof of this proposition is divided in several lemmas and is left 
to the end of this section. 

Let us first recall the co-area formula. Let 7r : X — > Y be a smooth 
surjective map between Riemannian manifolds X and Y. If almost 
every y G Y is a regular value of it, and <p : X — > K. is integrable, then 

<p(x)dX=[ [ ^-d^\y)dY. 

In particular, the co-area formula for the projection m : V — > H(d) and 
a function if : V — > E yields 

Proof of Theorem [H Recall from Proposition [H that (I) is defined by 

(I) c^ 3/2 /" /" /■ /*(/*, c* 



(27r)*voi(p(o+i)) y /6Ww) y C6P(c » +1) y t6[0il] ii/*ii 2 

•||n c (/)|| ||A(||cir d ')/(C)l|e- |l/l|2/2 rftdCrf/. 

Then, for < t < 1, by the co-area formula for the map $f : "H^) x 
P(C n+1 ) V, and Proposition&we obtain 

(1) (2vr)^vol(P(C"^)) A 7 (M)6V 1 



nn ( (/)ii mur'-mcn e . uf/2 d ^ r;) m dt 

If $ t (/,C) = (M), then /(C) = h(C)/t and = U c (h). From 

Lemma || we obtain that, for all t G (0,1), H t \B(h,ri) '■ B(h,r)) — > 
$ < T 1 (/i, rj) given by ( i-> (/i t (£),C)> is a parameterization of the fiber 
$ < r 1 (^ ) ? ?)- Moreover, since ( = H t (/, C) whenever i?t(C) = (/, C)> 
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applying the change of variable formula we conclude that 
(12) 



(I) 



CD 3 ! 2 



(2vr) 7V vol(P(C™+ 1 )) J 



t 



-2n-l 



MM; 

\h\\ 2 



NJ Wl (h,T])- 



n c (/i)|| ||A(||C||- dl )^(C)l|e- ||/lt(c)l|2/2 d(dVdt. 



'CeB(M) 

From the definition of h t {Q) in (ITT]) and the reproducing kernel property 
of the Weyl Hermitian product (jSJ), we obtain 



MOH 2 = IN 2 + 2 



1 - 1 



Re(h,A((C,Cr di (;O di )KO)+ 



1 -t\2 



|A«c,C>-*(-,C>*)^(C)ll ! 



and then, 
(13) 



i/ii 



|A(||Cir*)MC)ll : 



Therefore, from (IT2T) and (fTB"|) we obtain 



(14) (I) 
where 

e(M) = 



CD 3/2 
"(2^ 

1 



(h,»7)GV 



^^))wv,.H-iiA(iicinMoiir- 

•||A(||C||-*)/i(C)||-2;(||A(||C||-*)/i(C)|| 2 /2)dC, 
and X„(a) = e (i-*-> r 2n-i ^ 

Now, the proof of Theorem [I] follows applying the co-area formula 
for the projection 7Ti : V — > □ 

3.1. Proof of Proposition O The map h t : P(C n+1 ) — > Hu) given in 
(ITT]) is differentiable, and therefore if( is also differentiable. 

Lemma 5. Let (/, £) G "H^) x P(C n+1 ) such that $t is defined and let 

(M) = $ t (/,0- r/ien > 



det 



det(Jd f x + (D^(C)| c x)*.D^(C))| ( 



1/2 



NJ vl {h,rj) 



where [Id H{d) , -(Dh t (()\^)*) : ft (d) -> ft (<0 x T C P(C 
mapf^(f,-(Dh t (C)\ c ±)*f). 



•<n+l N 



is i/ie linear 
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det(a| 



Proof. In general, let E\, E 2 , and V be finite dimensional vector spaces 
with inner product. Assume that dim(V) = dim(£' 1 ), and let p : V — > 
Ei be an isomorphism. Let 7 : E 2 — > E\ and a : E\ x E 2 — > V be 
linear operators. Consider the following diagram: 

Eh x E 2 V 
.* 

(/d Sl -7*) .•' 

E\ E 2 E\ 

where (7, IoIe 2 ) ■ E 2 — > E\ x E 2 . 

Note that the image of the operator (IoIe 1 , —7*) : E\ E\ x E 2 is 
the orthogonal complement of (7, Id)(E 2 ) in EiX E 2 , therefore we get, 

_ |det (p ■ a ■ (Id El , ~7*))l 
l((7,M B2 )(^))^l J det ( IdjSi +7 . 7 *)|i/2. 1 det (p)| 

_ |det (p-a- (7d gl ,-7*))| 

= |det(Id S2 +7*- 7 )| 1 /2. |det(p)|' 

where the last equality follows by the Sylvester Theorem: if A and B 
are matrices of size n x m and m x n respectively, then 

(15) det(Id m + BA) = det(Id n , + AB). 

Now the proof follows taking E x = U( d) , E 2 = T ? P(C n+1 ), V = 
T(h,rj)V, with the associated real inner products, 7 = Dh t (C)\^-L, & = 
D$t(/,C) and p = D>ir 1 {h,ri). □ 

The derivative of h t at C e P(C n+1 ) in the direction ( e T f P(C n+1 ) 
is given by 

where K c , L c : T c P(C n+1 ) -»■ are given by 

(16) Jf c (C) = A(M^)-D/i(OC; 



1 ^4^f^ 1 



(C,C) f 
(CO' 



for all C G T ? P(C n+1 ). 

Lemma 6. T7ie adjoints operators K^*, L^* : "H^ — >■ T^P(C n+1 ), are 
gif en by 

(is) # c *(/) = (DHOi^r- A((c,c)- di+1 )/(o, 
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and 

(19) LCU) = (Df(0\^y-A(((,()- di+1 )h(0, 

for any f G H(d) ■ 

Proof. By the definition of adjoint, the definition of and the repro- 
ducing kernel property of the Weyl Hermitian product we get 

Re{K ( *(f),() = ||CI| 2 Re(/,A((C,C)- dl (-,C) dl ) ■ Dh{QC) 

= Re(/(C),A((C,C)^ +1 )-^(C)C) 

= Re(( J DMC)| c x)*-A((C,C)- d ' +1 )/(C),C)- 

Moreover, differentiating equation (jSJ) with respect to (, we obtain 
for L^* that 

Re(V(/),C) = ||C|| 2 Re(/, A((C,C)- d ^(-,C) dl - 1 (-,C))MO) 
= Re( J D/(C)C,A((C,C)" dl+1 )MO) 
= Re(( J D/(C)| c x)*-A((C,C)- d ' +1 )MC),C)- 

□ 

Lemma 7. One has, 

det(/d c x + (Dht(Ok±)* ■ Dht(C)\ c ±, 

i+(^) 2 iiA(^iicir di )Mc)ii : 

( (¥) Vmoi c o* ■ a (iicii-* +i ) 2 • Dh(ch ^ 



det 



Idfx i 

V 1 + l|A(V5iicii- d 0Mc)H 2 ; 

Proof. By direct computation we get 

K C * • JT C = WOlcO* • A «C C)^ !+1 ) ■ Dh{Q\ c i; 
K ( * ■ L c = L ( * ■ K c = 0. 

Note that, if / = L c (() for some C e T c P(C n+1 ), then, for all 9 G C n 
we get 

(Df(()\ c ,ye = (Re<0, a (rf 4 ||cir 2 ) MO)) C 

Hence, 



L C *L C 



MO 



Id cJ 
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Therefore we get: 



1 - t\2 



(Dh{C)\ c x)*-A 



-2d<+2 



)-Dh(0\^ 

l|-*U(C) 



Id, 



The proof follows. 
Lemma 8. One has 



□ 



det 



d(*i o $ 4 )(^(0, o ■ (idH (d) ,-(Dh t (o\ c ,y ^ 

= | det(/d c x + (D^(C)| c x)* ■ ^t(C)l^)! t 2n . 



Proof. First we find an expression for the term inside the determinant. 
For short, let 



^ = Dfa o $ t )(MC), • (Id W((0 , -(2?fc(C)| c x)*). 



One gets, 



(20) 

and 
(21) 



A (7ri0 $ t)(/)C) 



(/)=/-(l-f)A 







(CC> £ 



/(C), 



(7T l0 $ t )(/,C) 



-(1-t) 



(0 = 

A(^|i)-^/(CK + A 



(CO* 



/(C) 



Since h t (0(() = h(Q/t, and Df^OKOI^ = ^(C)lc^ from fl2Q 
and (1211) we get 



^(/) = / - (1 - 1) A 



■ 



(CO 1 

+ (1 - t) [ A (|^) ■ D K0k- ■ (Dh t (0k-Tf+ 



+ A 



(C0 £ 
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for all / G H(d\. That is, with the notation Kq and given in ([16 
and (JTTJ) , we get 

(22) 

V(/) = /-(i-t) 



A 



(CO 



1-t 



K C (K C * + L c *)f 



+ 



1 - t\2 



L C (V + V)/ 



for all f eU( d )- 



Note that if) = Id% (d) — L, for a certain operator L. Therefore 



det(V0 = det((Id-^ (d) — L)| Tin ^), where last determinant must be un- 
derstood as the determinant of the linear operator (Id^ (d) — L)| Ttti ^ : 
ImL — > ImL. 

The image of L is decomposed into two orthogonal subspaces, namely: 



A 



(■,0* 



a : a = (ai, . . . , a n ) T G C n > ; 



■<c,c>* 

i? f := {L ? H : w G T c P(C n+1 )}. 

Note that lmK ( = C c C ker L ? * and ImL c = R c C ker 
Consider the linear map 



r:C n ^C c , r(6)=AP 
Note that, r" 1 (a( 



<-,C>" 
,<C,C> £ 

IA 



■ 0' 



(C0 £ 



(CO* 

A(||CI|- 
II = IIA 



■A(||C||*)&, beC n . 

) ■ a. Since 



we conclude that r is a linear isometry between C n and C^. 
Let 



|A(V5||CI|-*)fc(C) 



Since 



A Jdi 



HO 



\w\ 



for all w G T^P(C n+1 ), we get that rj is a linear isometry between 
T ? P(C n+1 ) and i? f . 

Let lie ^ an d n^ c ^ be the orthogonal projections on and i?^ 
respectively 



24 



D. ARMENTANO AND M. SHUB 



Then | det(V0| is equal to the absolute value of the determinant of 

A 

C D y 

where A = t _1 oII c <.^|c' c ot, B = t^oUc^Ir^t], C = rj^oU^ipl^or 



and D = r] 1 o Il R i/j\ r or]. 



Straightforward computations shows that 



A = tldt 



t 



■aciicii 



-di+l- 



■ Dh({)U ■ (Dh(0\c 



A 



B 

C 



ia(>/*iicii-*)/*(c)iia(iicii- 



1 -t\2 



D= 1 



1 - t\2 



A(v / ^iicir di )Moii 2 ') m c x. 



Since D is invertible, we may write 

f A — BD _1 C B 
D 







B 




(o 


D 


hence det 1 


'A 


3 


Thus, 






|det(^)| = 


t 2n \ 


14 



/ 



1 -t\2 



A 



fc(C) 



2/i. 



det 



Idc« 



¥> a 



-* +1 )-^(C)lc- WC)lc-)*-A(||CI|- 



ii+lA 



V 



Observe that 



!+(¥) •||A(va-ncil-*)MOir 

* +1 ) a -D/i(0| fi = 

-* +1 ).i>/i(C)| fi )*.(A(||ci|-* +1 )- J DMC)l^ 



A 



Then, the proof follows from Lemma [7] and the Sylvester theorem ( TT5|) . 

□ 

Proo/ o/ Proposition The jacobian of ff t : P(C n+1 ) -»■ x P(C n+1 ) 
at C is given by 

det(Id^ + (Dh t (0\ c ±)*-Dht(Oy 3 



Then, the proof follows from Lemma [5| and Lemma [3 



□ 
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4. Proof of Theorem [2] and Theorem [3] 
4.1. Proof of Theorem El Recall that 

h (O = l|n c (/OII - ncWIWII* - n c (/0|| 2 /2), 

where X n ,(a) = /„* eCi-*" 2 )^- 2 ™- 1 dt. 

In order to prove Theorem\^we need some extra notation. Denote by 
S 2n+1 the unit sphere in C n+1 = R 2 ^ n+1 \ For all fceN, let a k : R k ^ 
K be the norm function er^rr) = ||x||. Given a measurable function 
let us denote by K R k((p) the expected value with respect 
to the standard Gaussian measure on M. k , i.e., 

Lemma 9. Let W be a complex vector space of dimension m with an 
Hermitian product (-, -)yy. Denote \\w\\ = \J (w, w)w Then, for every 
integrable function ip : R — y M. one has 



(2k) 



[ V(lkll)e- |H|2/2 ^ = E R2m (^o C r 2 , 

Jw&W 



Proof. If Wi, . . . ,w m is an orthonormal basis of (W, (•, -)w)i then, w±, 
. . ., w m , iwi, . . ., iw m is an orthonormal basis of the associated 2m- 
dimensional real vector space, namely W^, with the real inner product 
Re(-, -)w Define the linear map A : W H — > R 2m , by A(w k ) = e fc , 
and A{iwk) = e m+ k, where is the k-th element of the standard 
orthonormal basis of M 2m . A straightforward computation shows that 
A is a (real) linear isometry. Then, the lemma follows by the change 
of variable formula. □ 

Lemma 10. Let p > 1. Then, 

^H id) {Ph\\ P LV ) =E mN - n) (a p 2{N _ n) )E K2n (a p 2n l n p (a 2 2n /2)) . 

Proof. Given the canonical projection S 2n+1 — > P(C n+1 ), by the co-area 
formula we get 



( 23 ) IMIk = o uwrn+iK / d h {wfdw. 

2n vol(P((L ri+i )) J weS 2 n +i 

Recall that V w = {/ G T-L{d) '■ f( w ) — 0} and C w is the Her- 
mitian complement of V^, and that Ii w : H(d) Vw and H\c w are 
the orthogonal projections onto V w and C w respectively. Then we 
may write %(d) — Ku © C w . Denote by ip(a) = «I„(a 2 /2). Then 
9 h (w) = \\U w (h)\\ ^(\\U\ Cw (h)\\). Since \\h\\ 2 = \\U w (h)\\ 2 + ||n c J^)|| 2 , 
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by Fubini's theorem we get 

e -INI 2 /2 r e -||/ll 2 /2 r e -||9ll 2 /2 



(d) 



(2tt)^ Jfev w nn (^) N ~ n J 9 ec w KnnJ (^) n 



Since V w and C w are complex vector spaces with an Hermitian prod- 
uct of dimensions N — n and n respectively, the proof follows inter- 
changing in (1231) the integral sign with the sign of expectation, Lemma 
fj and the fact that vol(S 2n+1 ) = 27rvol(P(C n+1 )). □ 



Lemma 11. Let p > 1. With the above definitions, one has 

(i\ TfT , sfrrP \ — OP/2 r(A f -n+p/2) 

W % 2 <"-"H°2(JV-n)J - 1 T(N-n) ■ 

>p/2 r(r 
P T(n) 



ii) E K 2n (o"2 n • Z n p (cr 2 n /2)) < r< 'rtn/ 2 ^ > where the equality holds 



for p — 1 . 

Proof, (i): By definition, and integrating by polar coordinates, we have 

vol(^ 2(7V - n) -^ /"+ 00 



2(7V-«)-l^ V/2 d 

Then performing the change of variable w = p 2 /2 we obtain: 

E R2( iv-n)(a2(jv_ n )) = (271-)^-" J ' 

vol(5' 2(7V " n) " 1 )2 iV - n+p/2 " 1 



(27T) 



AT-n 



T(iV-n + p/2). 



The assertion follows from the fact that vol(5 fe ) = 2 r ^ fc+1 ^ 2 ^ . 
(ii) By definition of X n , we have 

(24) E K 2n « • In p (<m) = j't-^E^n (a p 2n • e* 1 -^/ 2 ) dt. 
Then, for fixed i G (0,1), 

| x ||P e ^-ii)ll-ll 2 /2 e -||x||V2 dx 
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Now, performing the change of variable u = (Jj — (p — l))p 2 /2, we get 
^2n e j - (2?r)n( _ ^ _ 1))n+p/2 1 {n + p/2) 



j.2n+p 



o P/2 r(n + p/2) ^ 



{ P -t 2 (p- l))«+P/2 r(n) 

Since t G (0, 1), and p > 1, we have the bound p — t 2 (p — 1) > 1, and 



hence 



(where the equality holds for p = 1). Then, the proof follows from 

□ 



Proof of Theorem The proof follows from Lemma [TO and Lemma 
ED □ 

4.2. Proof of Theorem [S 

Proposition 4. Let p, q > 1 s«c/i i/iai l/p+l/g = l. If 1 < q < 2, 

then we have 

~2P T(N -n + p/2) T( y n + p/2y 1/p 



(I) < C-D 3 ^ 2 D 



p T(N - n) T{n) 
T(N + 1) Y[n 2 + n-q) 2 2 "+ 2 3q 



V? 



_r(JV + l-g) r(n 2 + n) 4 - 2g 

Proof. Consider on V the following density measure: 

d Pv = (27iy N 'D~ 1 NJ 7ri (h,ri)e~ lM2/2 dV. 

By the co-area formula we get that this measure is a probability mea- 
sure on V. 

From (I14p . the definition of in ([7]), and Fubini's theorem, we obtain 

\2 



(i) = CD^j v ^-e(h, v )d P , 
<cdv*t> j^e^v 



voi(P(C"+ 1 )) 7 P(C „ +1) V7v IH 

The function 6h{(), as a function defined on V, is constant and equal to 
||n f (/z)|| ||/i-n c (/i)||X n (||/i - n c O)|| 2 /2) on the fiber of the projection 
tti : V^U [d) . 
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For p, q > such that 1/p + l/q = 1, Holder inequality on (V, dpy) 
yields 

Then, applying Fubini's theorem and taking the canonical projection 
7Ti : V — >■ "H(rf) we get by the co-area formula 

CDZI2 V U ( ^ /i 29 (M)" 



(I) < ^1/,+1/p E «( d) 



£)l/<?+l/p n M ^ M^ipg 

\v/Hv)=o 11 11 



vol(P(C» +1 ))i r(C „ + , 1 V m "« ( ''" (C) " ) ) 1/! '^- 



where, abusing notation, we denote 



^/h{rj)=0 """ ) ^ 1 Jhen (d) v /h( v )=0 

From the proof of Lemma [77j we get that (6h{() p ) is independent 
of C G P(C n+1 ). Then, 

®m (d) (e h (cy)=E n(d) (\\6 h \\i P ), 

and hence 
(25) 

(i) < cD*' 2r D l i p (e Hw ( ^M)) e ^)(II^II^) 1/p - 
V v/m^o / / 

In Bertram Shub (see also Beltran-Pardo [6]) it is proved that, for 
< a < 4, 
(26) 

£(M)^ < ^ r(jV+l) r(n 2 + n- Q /2) 2 Q + 2 ^ 3a/2 
- r(JV+l-a/2) T(n 2 + n) 4- a 




Now the proof follows from ( 1251) . ( |26i) and Theorem^ □ 

Proof of Theorem^ The Gamma function T(x), for x > 0, is logarith- 
mic convex (see Artin [JJ). Then, it is an easy exercise to check that 
T(x + 1/2) < y/xT(x), for all x > 0. 

Let p = 3 and q = 3/2. Then, applying this inequality to Proposition 
yields 

(J) < C'D 3/2 DN n y/N - n + 1/2 y/n + 1/2, 
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where C = C2 23 / 6 3 1/6 < 18C. Then, the theorem follows from the 
trivial bound N — n + 1/2 < N and n + 1/2 < (3/2)n. □ 

5. Numerical Experiments 

In this section we present some numerical experiments for n = 1 
and d = 7 that were performed by Carlos Beltran on the Altamira 
supercomputer at the Universidad de Cantabria. 

Recall from Theorem [I] that 

■iiA(iicir d 0Mc)iix ra (iiA(iicir d 0Mc)ii 2 /2)rfc, 

where X n (a) = J* e (i-*- 2 )<*t-2n-i dt. 

Table[l\ concerns a degree 7 polynomial h, chosen at random with the 
Bombieri-Weyl distribution. The condition numbers /j,(h,rj), Q(h,rj) 
and vol(B(h,rj)), at each root t] of h are computed. 

The data of the chosen random polynomial is given by: 



One gets \\h 



a-j = 


-0.152840 -iO. 757630 


a 6 = 


1.283080 + i0.357670 


a 5 = 


2.000560 + z3.302700 


CZ4 = 


13.004500 + z0.203300 


a 3 = 


-1.138140 + i7.094290 


a 2 = 


3.110090 + z2.618830 


ai = 


0.282940 + -z0.276260 


a = 


-0.316220 + ^0.036590, 


2.9631. 





Roots in C 


fi(h, ■) 


0(V) 


vol(5(M) 


3.260883 - il.658800 


1.712852 


0.4733570 


0.140509tt 


-2.357860 - il. 329208 


1.738380 


0.5502839 


0.138576tt 


-0.210068 + il.868947 


1.608231 


0.5049662 


0.144054vr 


0.227994 - i0.782004 


1.909433 


0.4914771 


0.125685tt 


-0.044701 + iO. 384342 


3.231554 


1.003594 


0.147277tt 


-0.308283 + i0.049618 


3.183603 


0.8892611 


0.152433tt 


0.213950 - i0.068700 


2.948318 


0.8426484 


0.151466tt 


Table 1. Degree 7 ranc 


om polynomial. 
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In i*%we H3 we have plotted, using GNU Octave, the basins B(h,rj) 
at each root t] of the chosen random polynomial h are plotted, in C 
and in the Riemann sphere,. 

Basins of atraction inthe plane Basins of attraction in the plane, zoom around and unit circle 




-0.8-0,6 -0.4 -0.2 0.2 0.4 0.6 0.8 -0.8 -0.6 -0,4 -0.2 0.2 0.4 0.6 0.8 



Figure 2. The basins B(h, 77) in C and in the Riemann 
sphere of the degree 7 random polynomial (GNU Oc- 
tave). 

In Table [H the same quantities are computed for the polynomial 
given by ao = —1, a± = . . . = = 0, aj = 1. In this case the roots 
are the 7th roots of unity, and it is not difficult to see that the actual 
values of fi(h,rj), Q(h,r]) and vo\(B(h, 77)) are constant at the roots of 
h by symmetry. This example illustrate the extent of accuracy of the 
computations. 

In this case we get \\h\\ = y/2. 

The errors for the root of unity case in the third column are of the 
order of 25%. But 25% does not seem enough to explain the variation 
in the computed quantities in the third column of the random example 
where the ratio of the max to min is greater than 2. So it is likely that 
they are not all equal. On the other hand, the ratios of the volumes 
of the basins in the fourth columns of the random and roots of unity 
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Roots in C 


MM 


0(V) 


vol(B(V)) 


-0.900969 + i0.433884 


3.023716 


0.7035899 


0.128982tt 


-0.900969 - i0.433884 


3.023716 


0.8354068 


0.153846tt 


-0.222521 + iO. 974928 


3.023716 


0.7405610 


0.135198tt 


-0.222521 - iO. 974928 


3.023716 


0.7549753 


0.1414147T 


1.000000 + i0.000000 


3.023716 


0.9128278 


0.156954vr 


0.623490 + i0.781831 


3.023716 


0.6800328 


0.135198tt 


0.623490 - i0.781831 


3.023716 


0.8122845 


0.148407vr 



Table 2. h(z , z\) — z\- z { 



Basins of atraction inthe plane Basins of attraction in the plane, zoom around and unit circle 




Figure 3. The basins B(h, rj) in C and in the Riemann 
sphere for h(zo, z\) = z\ — z 7 (GNU Octave). 

examples do seem of the same order of magnitude. So perhaps for n = 1 
they are all equal? Also, the graphics of the basins are very encouraging 
in the random case. There appear to be 7 connected regions with a 
root in each. So there is some hope that this is true in general. That 
is there may generically be a root in each connected component of the 
basins. This would be very interesting and would be very good start 
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on understanding the integrals. It would be good to have some more 
experiments and even better some theorems. 
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